LAGRANGE INTERPOLATION AT REAL PROJECTIONS OF LEJA 
SEQUENCES FOR THE UNIT DISK 



JEAN-PAUL CALVI AND PHUNG VAN MANH 

Abstract. We show that the Lebesgue constant of the real projection of Leja sequences 
for the unit disk grows like a polynomial. The main application is the first construction of 
explicit multivariate interpolation points in [—1, 1]"^ whose Lebesgue constant also grows 
like a polynomial. 



1. Introduction 

We pursue the work initiated in f4 | that aims to construct explicit and (or) easily com- 
putable sets of efficient points for multivariate Lagrange interpolation by using the process 
of intertwining (see below) certain univariate sequences of points. Here, the efficiency of 
the interpolation points is measured by the growth of their Lebesgue constant (the norms 
of the interpolation operator). Namely, we look for sets of points P„ C — for inter- 
polation by polynomials of total degree at most n — for which the Lebesgue constants 
A(P„) grows at most like a polynomial in n. We say that such points are good interpolation 
points in the sense that if A(P„) = 0{n") with a e N* then, in view of a classical result 
of Jackson, the Lagrange interpolation polynomials at P„ of any function with a + 1 con- 
tinuous (total) derivatives converge uniformly. This is detailed in the paper. In our former 
work IJl, multivariate interpolation points with good Lebesgue constants were constructed 
on the Cartesian product of many plane compact subsets bounded by sufficiently regular 
Jordan curves (including, of course, polydiscs) starting from Leja points for the unit disc 
(see below). Yet, from a practical point of view, especially if we have in mind applications 
to numerical analysis, the real case is more interesting. It is the purpose of this paper to 
exhibit explicit interpolation points in [— 1, l]'^ with a Lebesgue constant growing at most 
like a polynomial. As far as we know, this is the first general construction of such points. 
This will be done by suitably modifying the methods employed in Q. Actually, the uni- 
dimensional points will be taken as the projections on the real axis of the points of a Leja 
sequence. We shall first show how to describe (and compute) these points and, in partic- 
ular, prove that they are Chebyshev-Lobatto points (of increasing degree) arranged in a 
certain manner. We shall then study their Lebesgue constant to prove that it grows at most 
like rr' log « where « is the degree of interpolation. The passage to the multivariate case is 
identical to that shown in f?) and will not be detailed. 

Notation. We refer to [4| for basic definitions on Lagrange interpolation theory. Let 
us just indicate that, given a finite set A, we write wa ■— YlaeAi' ^ '^)- The fundamental 
Lagrange interpolation polynomial (FLIP) for a e A is denoted by £{A,a; •). We have 

t[A,a--)= n {■-b)/{a-b) = WA/{w'A{a){--a)). 

beA.b^a 

The Lagrange interpolation polynomial of / is L[A;/] = Y.aeAf{'^)^{^i'^'^ ') the norm 
of L[A;-] as an operator on C{K) (where K is a compact subset containing A) is the 
Lebesgue constant A(A) = A{A,K). It is known that A(A) — max,^;^:^^^^ |^(A,fl;z)|. 
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We shall denote by M, M', etc constants independent of the relevant parameters. Occur- 
rences of the same letter in different places do not necessarily refer to the same constant. 

2. Leja sequences and their projections on the real axis 

2.1. Leja sequences for the unit disk. We briefly recall the definition and the structure 
of a Leja sequence for the unit disk D = { |z| < 1 } C C. A A;-tuple -Ejt = (eo, ■ • • , Sk- 1) & 
with eo = 1 is a k-Leja section for D if, for j — l,...,k~l, the ( j + 1 ) -st entry e j maximizes 
the product of the distances to the j previous points, that is 

J-l .7-1 

n |ey-em| = max f| \z-e„\, J = l,...,k-l. 

The maximum principle implies that the points e, actually lie on the unit circle dD. A 
sequence E = {e/t : k E N) for which :~ (eo, • • • ,^k-i) is a A:-Leja section for every 
^ e N* is called a Leja sequence (for D). Of course, the points of a Leja sequence are 
pairwise distinct. 

The structure of a Leja sequence is studied in H] where one can find the following 
result. 

Theorem 2.1 (Bialas-Ciez and Calvi). A Leja sequence is characterized by the following 
two properties. 

(1) The underlying set of a 2" -Leja section for D is formed of the 2"-th roots of unity. 

(2) If E 211+ 1 is a 2"^^-Leja section then there exist a 2"-root p of — 1 and a 2"-Leja 
section £'2''' such that E-,1,+1 = {Ej" , pE^,^). 

Repeated applications of the above rule show that if t/ = 2"" + 2"' H h 2"'' with no > 

«!>•••> > then 

Ed = (Ero , Po^^Ly-o) = (Ero , Po^y'l ' PiPoEd'-2"Q-ri ) ^2.1) 
= • • • = {E2"o , PqE^I] , piPo£'2"2 > • • • , Pr-i ■ ■■PiPoE2nl), (2.2) 

where each E^Jj consists of a complete set of 2"' -roots of unity, arranged in a certain order 
(actually, a 2"-* -Leja section), and pj is a 2"j-th root of — L 

2.2. Projections of Leja sequences. We are interested in polynomial interpolation at the 
projections on the real axis of the points of a Leja sequence. We eliminate repeated values 
and this somewhat complicates the description of the resulting sequence. We use to 
denote the real part of a complex number (or sequence). 

Definition 2.2. A sequence X (in [—1,1]) is said to be a Si-Leja sequence if there exists a 
Leja sequence E = (ei^ : G N) such that X is obtained by eliminating repeated points in 
3i{ek :k<EN). We write X =X{E). 

In particular X(£') is a subsequence of3i{eii : k eN). Since for every « e N, the under- 
lying set of a 2"+^-Leja section is a complete set of 2"+'-st roots of unity (Theorem 



2.1 



the corresponding real parts form the set of Chebyshev-Lobatto (or Gauss-Lobatto) 
points of degree 2", 

if2"-{cos(;;r/2") : j = 0,...,2"}. 

These points are the extremal points of the usual Chebyshev polynomial (of degree 2") and 
are sometimes referred to as the "Chebyshev extremal points". 
For future reference, we state this observation as a lemma. 



Lemma 2.3. Let X be a 3i-Leja sequence. For every n G N, the underlying set of X2"+i — 
(xo, . . . ,X2") is the set of Chebyshev-Lobatto points ^2"- 
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Theorem 2.4 below gives two descriptions of 9{-Leja sequences. The first one is partic- 
ularly adapted to the computations of 9{-Leja sequences when one is given Leja sequences. 
Examples of easily computable (and explicit) Leja sequences can be found in |4 Lemma 
2]. In Figure [T] (I), we show the first 16 points of a Leja sequence E and the first 9 points 
of the corresponding 5R-Leja sequence X{E). (The Leja sequence we use is given by the 
rule £2 = (1.^1) and/iji'+i = (£'2'',exp(/;r/2")£'2")-) 

The concatenation of tuples is denoted by A, 

(xi, . . . ,x,„) A (yi , . . . ,y„) := (xi , . . . ,x„,,yi , . . . ,y„). 

For every sequence of complex numbers 5 = (ijt : G N) we define S{j :k): = {sj,Sj+i ,Sk) 
As before, 5i::= 5(0: 

Theorem 2.4. A sequence X = [xj^ : k ^H) is a 3i-Leja sequence if and only if there exists 
a Leja sequence E = (ei^ : k CzN) such that 

X = {1,-1)A /\3{(e{2' -.ij + 2j-^ -l)y (2.3) 

Equivalently x^ — 9{(e||,(^)), A: G N, with 0(0) = 0, 0(1) = 1 and 

<l>(k) = < \ /ni , k>2, (2.4) 

> |2Li°g2WJ k^2" - 

where [-J is used for the ordinary floor function. 

Proof Let X = X{E) with £ = (e, : 5 G N). We prove that if V <k< 2'i +2'-^ then 
5R(e<;) does not appear in '^{Ej^) and therefore provides a new point forX. To do that, since 
Ck itself does not belong to Ek, it suffices to check that 'e^ is not a point of Ek, equivalently 
S/t 7^ ^1 < i < fc — 1. If 5 < 2-' then is a 2''-th root o f un ity whereas ei, is not. On the 
other hand, if 2^ < < — 1 then, in view of Theorem 2.1 = pa and Cs = pb where 
p is a 2''-th root of —1 and both a and b are 2''^'-st roots of unity. The relation e^ =e^ 
yields p /'p = b/a. The argument of the first number is of the form 2(2Z + l);r/2'' and the 
argument of the second on e is 2 f;r/2-^~' (with / , f G Z). Equality is therefore impossible. 



Now, in view of Lemma 2.3 from £2;+! we obtain 2^ + 1 points forX, namely the points 
in arranged in a certain way. Yet, the 2^ + 1 first points of X are already given by 
(2-'^' + 1 points) together with, according to the first part of this proof, the 2^^^ points 
^{ek), 2j <k< 2^+2'-i. This implies that ifV+V-^ <k< 2^+1 then 3i{ek) is not a 
new point for X. This achieves the proof of ( |2.3| l. 

To prove ( |2.4| i we observe that, in view of p.3| l, we have 

9^(e2.+,)=^2*-'+,+i. 0</<2*-'-l. 
Hence (2*^^ ' + / + 1 ) = 2'-' + / and the expression for (j) easily follows. □ 

CoroUary 2.5 (to the proof). IfX = X{E) then X{2" + 1 : 2"+') = 3{{E{2"+^ ; 2"+' +2" - 
!))■ 



Decompositions ( |2.2| l and ( |2.3| l are fundamental to this work. In particular, the binary 
expansion of k will be used in the study of the tuple X{Q : k). 

Note that, of course, the decomposition would be different if we projected the Leja 
points on another segment, say on [—e'^,e'^]. The distribution of the projected points in 
general depends on arithmetic properties of 9. We shall not discuss the general case in this 
paper 

Finally, let us point out that our 9{-Leja sequences are not Leja sequences for the in- 
terval. It can be shown that they are pseudo Leja sequences (see 1 1 1 for the definition of 
pseudo Leja sequences). There is no known expression for Leja points for the interval. For 
that matter, such an expression is very unlikely to exist. 
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(I) First 9 points of a 9?-Leja sequence. (II) 45 interpolation points obtained as the intertwin- 

ing of the points in (I) with themselves. 

Figure 1 . Exemples of points from a 9{-Leja sequence and their intertwining. 

3. LEBESGUE constants of 5R-LEJA SEQUENCES 

3.1. The upper bound and its consequences. Recall that given a set A of n + 1 inter- 
polation points in [—1,1] and / e C'*'([— 1, 1]), the Lebesgue inequality together with the 
Jackson theorem [7 Theorem 1.5] yield the well known estimate 

max ||/-L[A;/]|| =M(1+A(A)) co{f^'\2/n)/n' 

where a)(/^*' , •) denotes the modulus of continuity of /''^^ and M does not depend either on 
A or n. The following theorem implies in particular that interpolation polynomials at the 
points of any 9^-Leja sequence converge uniformly on [—1, 1] to the interpolated function 
as soon as it belongs to C'*([— 1, 1]). A weaker consequence is that the discrete measure 
l^d ■~ 3TT ^f=o[-*^'] associated to the SR-Leja sequence X = {x^ : s ^ N) weakly converges 
to the 'arcsin' distribution on [—1,1] which is the equilibrium measure of the interval. Here 
[xi] denotes the Dirac measure at jc, . 

Theorem 3.1. LetX be a 3i-Leja sequence. The Lebesgue constant A(Xii) for the interpo- 
lation points xq,. . . ,JCi;_i satisfies the following estimate 

A{Xk)^0{khogk), ;t^oo. 

The construction of good multivariate interpolation points is derived as follows. We 
start from 9{-Leja sequences X'^.') — {x[^^ :kEN),j = l,...,N. These A^ sequences need 
not be distinct. We define P<. c [— 1 , l]'^ as 

P^ = {x„-(4\...,xL');^):f a,<4- 

It is known 13 1 that this is a valid set for interpolation by A^-variables polynomials of degree 
at most k. Actually, Pi^ is the underlying set of the intertwining of the univariate tuples 
X'^j\0 : k), j = I,. . . ,N. We refer to g) for details on this definition and to |3| for a 
general discussion of the intertwining process. Let us just emphasize that, in order to 
provide good points, the method requires to use sequences of interpolation points (we add 
one point when we go from degree kto k+1) rather than arrays (the classical case : all the 
points change when we change degree). The reason for this is explained in |4|. 

In Figure[T](II), we show the points of a set constructed with the first 9 points of the 
5R-Leja sequence in (I). 

Theorem 3.2. The Lebesgue constant A(P<;) grows at most like a polynomial in k as k ^ °°. 
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Proof. The proof of |4, Theorem 16] works as well in this case. □ 
Just as in the univariate case, the multivariate versions of the Lebesgue inequality and 



the Jackson theorem together with Theorem 3.2 imply that Lagrange interpolants of suf- 
ficiently smooth functions converge uniformly. Examining the terms in the proof of ID 
Theorem 16], we find 



A(P,)-o(fe(^'+'i'^-'^)/2log'^fe), k 



which gives a more precise idea of the required level of smoothness. This bound however 
is certainly pessimistic. 



3.2. Outline of the proof of Theorem 3.1 We take advantage of the sti'ucture of the 



points of a 9{-Leja sequence. The first step is a simple algebraic observation. We use the 
notation recalled at the end of the introduction. 



Leiimia 3.3. Let N 

We have 

i{N,a;-) = 



NqU ■ ■ ■ UNs-i where the Niform a partition of the finite set N C 



WN\Ni{a) 



i{Ni,a\-), a<ENi, / = 0, 



Consequently, 



s-l 

MN) < y max 
^^xeK.ueNi 

!=0 ' 



(3.1) 



(3.2) 



WN\N,{a) 

where the Lebesgue constants are computed with respect to a compact set K containing N. 

Proof. We readily check that the polynomial on the right hand side of ( |3.1| l satisfies the 
defining properties of £{N,a; ■). The estimate for the Lebesgue constant A(A^) follows from 
the definition together with the formula for the FLIPs. □ 

Given a 9^-Leja sequence X ^ X{E), to estimate A{Xk), 2" + 1 < k < 2"+', we shall 
first apply the lemma with a partition of X/^ into two subsets, namely A = X{Q : 2") = 
and B = X{2" + I : k — 1). Here and below, when there is no risk of misunderstanding, 
we confuse a tuple with its underlying set. In other words, when we write T — Y with 
T a tuple and Y a set, we mean that Y is formed of the entries of T. Our choice for A 
and B, of course, is motivated by the fact that the Chebyshev-Lobatto points are excellent 
interpolation points for which there is a large amount of information available, see below. 
With this choice. Lemma [33] gives 



A(Xi) < A(A) max 

xg[-i,i1.«ga 



wb{x) 



wsia) 



-A(B) max 

xe\-i.i].heB 



wa{x) 



WAib) 



(3.3) 



The factors depending on the first subset A in p.3| l will be easily estimated. The more 
difficult part will be to estimate wg and Ag. To do that, we shall use a partition of B 
and still have recourse to Lemma 3.3 in its most general form. We point out however 



that our method is unlikely, it seems, to give the best estimates. Intuitively, we think that 
sharp estimates cannot be obtained by separating the interpolation points into two or more 
groups. 

It is not difficult to see that the Lebesgue constant A{Xii) cannot grow slower than k. 



This is explained below in subsection 3.4 



3.3. Interpolation at Chebyshev-Lobatto points. We collect a few results on Chebyshev- 
Lobatto points. First, since they are the extremal points of the Chebyshev polynomials, we 
have 

w^^{x)^{x^~l)d-'T;{x), xeM, 

where denotes the monic Chebyshev polynomial of degree d. From 2'^^^Td{cos9) = 
cosdO we readily find 

w^^(cos0) = -2'"''sin0sin£/0, G M. (3.4) 
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A classical result of Ehlich and Zeller 15i|2J ensures that 

Ai^^)=Oilogd), d^^. (3.5) 

Lemma 3.4. Let X = X{E) be a 3i-Leja sequence. If2" + 1 <k< 2"+', A=X{Q: 2"), 
B=X{2" + 1 : 1) and K= [-1,1] then 

max ^<l/sin(;r/2"+'). 

xeK,beB \wA{b)\ 

Proof. If X = cos t and b = cos j3 then, since A = , in view of p.4| l, we have 
\wA{x)/wA{b)\ = |sinfsin(2"f)|/|sini3sin(2")3)| < 1/| sinj3 sin(2")3)|. 
Hence, it suffices to check that 

|sin/3sin(2"j3)| > sin(;r/2"+i). (3.6) 
Since b — xj with 2" + 1 < j <k—l < 2"+^, equation p.4| i gives 

x,=9^(e2'.+i-i)=5R(e2«+i+„_i), u^j-2", 1<m<2"-1. 



Theorem 



2.1 



says that e2"+i+M-i ~ where p is a 2"+'-st root of —1 and z is a 2"-th root 
of 1 . This means that the angle j3 that gives b = xj is of the form 

J3 = (2T+l);r/2"+i+2T';r/2" withT,T'eZ. 

It follows that I sin 2" i3 1 ~ 1 and I sin PI > sin(;r/2"+i). This gives inequality p.6| ) and 
concludes the proof of the lemma. □ 

3.4. A lower bound. We now show that when we remove one point from the Lebesgue 
constant grows significantly faster. Write a; = cos{in/d). Suppose that Aj := ^/ \ {aj} 
with I < j < d — I {so that we remove a point different from 1 and —1). We compute the 
values of the FLIPs for Aj at the missing point aj . From w_^j ( • ) — waj{-){- — aj) we get 

WAjiaj) = w'y,^{aj) and w'j^^{ai) ^w'A.{ai){ai~aj), i^j. 

Hence, 

Yet, as easily follows from p.4| l, 

w^^(a,) = ±2^^'^'c/, / = 1, . . . 1. 

Hence |£(A,-,fl,;flj)| = 1 for d — 2 values of /, namely for / = 1, . . . ,ii — 1, i ^ j. Conse- 
quently, 

d-l 

A(Ay)> \e{Af,ai-aj)\=d~2. (3.7) 

Here is the consequence about our 9{-leja sequences. If X is a 5R-Leja sequence, then 
X{Q :2" — 1) is formed of all the Chebyshev-Lobatto points of degree 2" with only one 
missing and this missing point is different from 1 and —1 (as soon as « > 2) which are 
the first two points of the sequence. Hence, according to p.7[ ), A{X{0 : 2" — 1)) > 2" — 2 
which shows that the Lebesgue constant A{Xii) cannot grows slower than k. 
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3.5. Interpolation at modified Cliebysliev points. We intioduce other sets of interpo- 
lation points that will naturally come into play when dealing with B = X{2" + 1 : k — I). 
When cosjS is not an extremal point of Tj, that is cos/3 ^ then the equation Td{x) = 
r^(cos/3) has d roots in [—1. 1]. The set of these roots — which we call modified Cheby- 

shev points — will by denoted by -3^}^^ Since the change of variable x = cosf transforms 
the equation in cos dt = cosdjS, we have 

£r}P^ = {cospj : pj := p +2j7t/d, j ^0, . . . ,d ~ I}. 
The following result is probably known but we are unable to provide references. 

Lemma 3.5. We have A{.^j^^) = 0(logt//| sint/j3 1) as d ^ °° and the constant involved 
in O does not depend on ji. Equivalently, there exist M such that. 



Proof. First, since 



A(J^/0 <log('^+l)/|sin£/j3|), d>l. 



w {x) = Td (x) - Td (cos /3 ) , 



see the introduction for the notation w (^), we have 

Td {x) — Td {cos j5) sin j5j {cos dt — cos dj5) 



e{,7}'^\cos j5j;x) 



rj(cos/3j) {x — cospj) d sind P {cost — cos Pj) 



■ cosf. 



A use of the sum-to-product identity for cosines now yields 



^(^/\cosj3;;^) = 



sinPj sin(t/(f +j3)/2)sin(J(/-j3)/2) 
dsindP sin((f +j3/)/2)sin((f-)3y)/2) ' 



(3.8) 



Yet, we also have 

sinj3;=sin((f + j3,)/2)cos((f-j3,-)/2)-sin((f-j3,-)/2)cos((f + j3,)/2). 
Using this in p.8| ), we obtain after simplification, 

cos((f - l5j)/2) sm{d{t + 15) /2) sm{d{t - j3)/2) 



^(^/\cosj3^;x) = ^ 



dsin dp 



sm{{t-pj)/2) 
cos{{t+Pj)/2)sin{d{t + p)/2)sm{d{t-P)/2) \ 
sm{{t + pj)/2) /• 



(3.9) 



It follows that 



\£{£rjl^\cospj;x) \ < 



1 



d\smdp\ 



sm{d{t-p)/2) 



sin((f-i3,)/2) 



sm{d{t + p)/2) 



sin((f + /3;)/2) 



cosf. 



Hence, 



sin£/j3| 



maxF(f-j3)+maxF(f +j3) > = 



\smdp\ teR 



maxF(f), 



where 



F{t) 



d-l 



j=0 



sm{dt/2) 
sm{{t-2jn/d)l2) 



But max,gRF(f ) is exactly the Lebesgue constant for the li-th roots of unity which is known 
tobe C>(logt/), see Hg). □ 
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4. Proof of Theorem I3.1I 

4.1. Further reduction. We use Lemma [34] and the classical estimate ( |3.5| l of Ehlich and 
Zeller in p.3|l to obtain the following lemma. 



Lemma 4.1. LetX be a 3i-Leja sequence and let2" + l<k< 2"+K If A = X(0 : 2") and 

B^X{2" + 1 :k-l) then 

A(X,)<Mlog2« max ^ + ^^^^, (4.1) 
A'G[-i,i],aGA IwB^fljl sm(;r/2"+^) 

where M does not depend on k. 

The remaining required information is collected in the following two theorems. 

Theorem 4.2. LetX be a Si-Leja sequence and let!" + I <k<2"+\ If A ^ X {Q : 2") and 

B=X[2" + \ -.k-l) then 

max ^<22"+2. 

Theorem 4.3. Let X be a 3i-Leja sequence and let T^\<k< 2"+'. IfB = X{2" + \: 
k—\) then 

A(B) <M'22"log2", 
where the constant M' does not depend on k. 



End of proof of Theorem \3.1\ When ^ — 1 is a power of 2, the points of Xj^ form a complete 
set of Chebyshev-Lobatto points and the bound is imp lied by Ehlic h and Zeller's estimate 
p.5|l. We assume 2" + 1 < k < 2"+' . Using Theorems 



4.2 



and 



4.3 



in (14. IL we obtain 



AiXk) < M22«+2log2" +M' . f (4.2) 

sm(;r/2''+') 

Since n = [logzWJ (or [logsWJ - 1 in the case k = 2"+^) and l/sin(7r/2"+i) = 0(2"), 
this readily gives the existence of a constant M" (independent of k) such that 

A{Xk) < M"k^ \ogk, for k lai-ge enough. 

Observe that the highest power (that is, k^) comes from the second term in (|4.2|l. □ 



4.2. A trigonometric inequality. The proofs of the two remaining steps rest on an ele- 
mentary inequality that we present in this subsection. As in [4|, the key observation is 

|sina| > |sin2"a|/2", n e N, a e M. (4.3) 

Lemma 4.4. Let r > I and let no > n\ > ■ ■ ■ > nr > be a finite decreasing sequence of 
natural numbers. If2"](pj — % [27r] {i.e. 2"j(pj = % mod 2%), y = 0, . . . , r — 1, then 

r-l 

]~[|sin2"'+i-^(9-9o (Pj)\>{l/2"o-"')\cos2"o-^(p\, (p eR. (4.4) 

j=o 

Proof. The proof is by induction. To treat the case r = 1, we prove that 

\sm2"'-\(p-(po)\ > (l/2"o-'")|cos2"o->|. 
Using ( |4.3[ ) with a = 2"'^' — (po) and n = n^ — n\ we obtain 

|sin2"'-'((p-(po)| > (l/2"«-«')|sin2«o-i((p-(po)|. 
But, since 2"''(po = ?r [2n], \ sin2"0"'((p — (po)| = |cos2"0"'^| and the claim follows. 
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We now assume that the inequality is true for r = k and prove it for r = fc + 1. The 
induction hypothesis applied to ^ — ^ instead of <p yields 

k 

n|sin2'V+i-i(((jr)-,po)-<Pi (P;)l>(l/2"'""'^+')|cos2"i-i((p-(po)|, (p e M. 

;=i 

(4.5) 

multiplying by the term corresponding to j = 0, we obtain 

n|sin2"'+'-i(((p-(po)-(Pi (P/)| 

^ ^^^i^ I '^"2"' - ^ ((p - (po) cos 2"' - ' ((p - (po) I 

= ^;;rWTrl^i"2"'((p-(po)|, (peM. (4.6) 

Another use of ( |4.3| l with n = «o — n i — 1 shows that 

2«.-L. + i Nin2"' ((p - (po) I > I sin 2«o- ' ((P - (Po) I , <P G M. 

The sine on the right hand side is shown to be | cos 2""^'^] as in the case r — I. □ 



4.3. Proof of theorem |4.2[ Let X = X{E) and B = X{2" + \:k-\) with 2" + 1 < /t < 
2"+i. We write 

=2"+2"i +...+2"' withn-1 >«!>..• >nr>0, (4.7) 
and, to simplify the notation, 

no = n + l, (4.8) 
c/,. = 2"o + .-. + 2"', / = 0,...,r. (4.9) 

Then, in view of Corollary [23] we have 

r-l 

B = X{2" + 1 : 2" + (t/i - da)) A /\ X(2" + (t/; - do) + 1:2" + - do)) 

i=i 

r-I 

= /\9?(£(fl',:fl',+i-l)). 

!=0 

Now using the structure properties of a Leja sequence, see Theorem |2.1| and \2.2\ , we see 
that the points of X(2" + 1) are certain modified Chebyshev points (see subsection 



33|. Indeed, for / e {0, . . . , r - 1 }, 

E{di : di+i - 1) - p, • • -pipo^;,;;/, with = -i, (4.io) 

3i{E{di : di+i - 1)) = ^^tt'"^^'^ with j3, = argp,- = (2f,- + l)n/2"\ tj e Z. (4.11) 
This implies the following relation for the polynomial wb, 

^Bix) = Y[ 1 ^2";+i (-*^) - 72"'+i (cos(j3o H h A)) f ■ 

i=0 ^ 

It follows that for fl = cos 9 e A = X(0 : 2") 

max ^ . max'ri I cos(2"---0 - cos(2»..Kfi) + ■ ■ ■ + (^.,2) 
.vG[-i,i] |wB(fl)| rGM |cos(2"'+i9)-cos(2"'+i()3o + ---+A))l 

r-I 

< 2'-n 1/1 cos(2«'+'(p) -cos(2"'+' (/3o + - • • + (4.13) 
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Now, a use of the sum-to-product formula for cosines together with two applications of 
Lemma 4.4 (first with (jO; — jS,, then with 9, = — j3,) enable us to bound the denominator in 
( |4.13| l and arrive to 



max ; -— T < 



^G[-i,i] kB(fl)| ~ cos2(2"o-i^)' 

It remains to recall that «o = n + 1 so that 22(''o^"'') < 2^"+^ and observe that, since A = ^2", 
2"o-i(p = 2"(p = \%\ so that cos2(2"0"'^) = 1. This concludes the proof of Theorem 



4.2 



4.4. Proof of theorem |4Jj We still use the fact that, for X = X(£') and 2" + 1 < < 
2"+', we have B = Apo's,-, where the underlying set of B,- is with /3; = (2f,- + 

l)7r/2«^f^'eZ. 



max 



Using first Lemma 3.3 (with A^, = B,) and then Lemma 3.5 to bound A(By), we obtain 

^MbWI 102f2"'+l + l) 

A(B)<My max ^ ^ . -i^^^^^ „ (4.14) 

'- jeo-ehi.iJ.^eB, |wB\B/fl)| |sin(2«.+ i(A) + ---+j3,))| 

Now, just as in ( |4.13[ ) (we just remove one factor), for a — cos Qj E Bj, we have 

\^B\B,ix)\ ^^^^ |cos(2"-+'0-cos(2"--+'(A) + ---+/3,))| 

f-Ti] |wB\B/fl)| ?gr'',J|^^.|cos(2«.-+i0,-)-cos(2«.-+i(/3o + ---+A))I 

<2'-i n l/|cos(2"'+'0,)-cos(2"'+'(j3o + ---+A))l- (4.16) 

Again the sum-to-product formula for cosines transforms the right-hand side in a product 
of sines and it follows that the j-term in the right-hand side of ( |4.14| i is bounded by the 
maximum when a = cos Oj runs over Bj of 

log (2"^+ ■ + 1 ) n I sin" ' 2«'>' - ' (0; - j3o A) ifl I ^i"" ' 2"'+' - ' ( 0,- + j3o + • • • + A) | 

1=0. i/j 1=0 

(4.17) 

Here we used the fact that 

I sin2"^+i -1 (0, + Ao + • • • + A,) I = I sin2"^+i (Ao + • • • + A;) I (4. 18) 

which enabled us to insert the isolated sine in ( |4.14| i into the second product of ( |4.17| l. To 
prove ( |4.18| ), we observe that, since a ~ cos Oj E Bj, we have 

2"^+' 0; = 2"^+i (Ao + • • • + A,,) [2^] . (4. 19) 

We now estimate independently both products in ( |4. 17| i. The same bound is valid for 
every a E Bj and therefore provides an upper bound for the maximum over Bj as required. 

I) We start with the first product. In view of ( |4.19| l, since n,+i > rij+i whenever / < j, 
we have 

2"'+' (0, - Ao A) = 2"'+' (A,+ i + • • • + A>) [2?r], < / < 

On the other hand, since 2"'A4 — 7i [2;r], we also have 

2"'+i(A,+i + ---+Ai) = ?r[2;r], 0</<;. 
Thus the absolute value of the first j sines equals 1 and we just need to estimate 

'fj |sin2«'+i-i(0,-Ao A)l- (4.20) 



LAGRANGE INTERPOLATION AT REAL PROJECTIONS OF LEJA SEQUENCES FOR THE UNIT DISK 



II 



To do that, we apply Lemma 4.4 with (p = Oj — j5o — jij. We obtain the lower bound 

( 1 /2'0+i ) I COS 2"J+' -\ej-j3o j3y) I . Yet in view of ( |4l9l ) this cosine equals ± 1 

and we obtain 

r-I 

Yl \sm-^2"-+'-^i9j-l5o j3,)| <2"j+'-"\ (4.21) 

Note that, in the case j = r—l, the whole product equals 1 which obviously implies the 
inequality. The inequality is likewise satisfied in the case r = 1 (for which the product is 
empty). 

II) We now turn to the second product in ( |4.17[ ). Using again Lemma 4.4 with (p = 
9j + /3o as in ( |4.6| l, we get the upper bound 

2"i-«,-+i/|sin2"'(0^'+j3o)|. 

However, since for every s, )3,. = (2f,. + l);r/2"-' and 0/ = )3o H h +2qjn/2"j+^ with 

ts , qj E we have 



2"' (Oj + M^ 2"! (2j3o + j3i + ■ ■ ■ + /3, + 2^,-;r/2"^+i ) 

= 2'"+% + p;r=^g±^+p;r withpeZ. 



Note that since no = « + 1 and n>n\ we have riQ — ni — l > 0. This shows that | sin2"i (0^- 
j3o)| > sin(;r/2"o-"i-') > 2/2"o-"i-^ = i/2"o-«i-2. We have therefore proved 



r-l 



]^|sin-i2"'+i-'(0j+j3o + ---+AO| <2"'""''+'-2"o-"'-2^2"o-"'-i <2". (4.22) 

III) It remains to insert ( |4.21| l and ( |4.22| l in ( |4.14| l with the aid of ( |4.17| l. Indeed, we 
obtain 

r-l 

A{B) < M ^ 2"J+'-"'- ■ 2" ■ log(2"'"+i + 1) (4.23) 

.7=0 

r-l r-I 

< M2" £ 2"^+^ log(2"j+' + 1) < M2"log(2" + 1) ^ 2"^+' = O(22"log2"). (4.24) 

;=o ;=0 

This achieves the proof of Theorem |4. 3 [ 
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